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Recently NA60 collaboration has extracted the inverse slope parameters from the dimuon spectra 
originating from In+In collisions at ^/snn = 17.3 GeV for various invariant mass region. They 
have observed that the inverse slope parameter as a function of invariant mass of the lepton pair 
drops beyond the p-peak. In the present work we study the slope parameters for the low mass 
region (LMR) and show that the experimentally observed behavior can be explained by assuming 
the formation of a partonic phase which reverts to hadronic phase through a very weak first order 
(or continuous) phase transition at a temperature T c ~ 175 MeV. A strong first order transition 
fails to explain the data. It is observed that a scenario without the formation of a partonic phase 
does not reproduce the non-monotonic behaviour of the inverse slope parameter non typical of radial 
flow. It is argued that the variation of T c g with invariant mass M may be used to constrain the 
transition temperature. 
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Low mass lepton pairs are considered as useful 
tools to probe the thermodynamic state of the mat- 
ter presumably formed in nuclear collisions at rel- 
ativistic energies (T), H, 0, 0, [E @ ■ Unlike hadrons 
which bring the information of the state of the mat- 
ter when they cease to interact - the lepton pairs 
probe the entire stage of the evolution history - from 
the initial formation time to the final freeze-out time 
of the fireball. Recently high quality data on dimuon 
spectra from NA60 collaboration [7| has been made 
available for both the kinematical variables - the 
transverse momentum (pt) and the invariant mass 
(M) of the pairs. They observed a non-monotonic 
behaviour of the inverse slope parameter, T c g, ex- 
tracted from the transverse mass spectra of the pairs 
- as a function of invariant mass and indicated that 
such a trend may be originated due to the radiation 
of lepton pairs from a partonic phase formed initially 
in the collisions. 

The probability that a muon pair of invariant 
mass, M and transverse momentum, px will emit 
from a thermal system at temperature T is deter- 
mined by the Boltzmann factor ~ %/ M2 +Pr/ T . For 
a dynamically evolving system like the one produced 
after nuclear collisions at ultra-relativistic energies - 
the temperature decreases with time because a part 
of the thermal energy is spent to allow the collec- 
tive motion of the system. Consequently the Boltz- 
mann factor is modified to ~ e -^/ M2 +PT/ T >=« ; where 
Toff ~ ?th + MuJ - here the first term represent ther- 
mal part and the second term stands for the flow part 
for flow velocity v r . It is expected that the large M 
thermal pairs will be emitted from early time when 
temperature is large and flow velocity is small and 
the small M pairs will originate from the late stage 
of the evolution when the temperature is low but 



flow velocity is large. Therefore, the variation of the 
T e ff with M may be treated as a chronometer of the 
heavy ion collisions. 

In the present work, therefore, we focus on the 
transverse mass distribution of the lepton pairs 
formed in In + In collisions at 158 A GeV beam en- 
ergy at CERN-SPS and subsequently extract T G fr- 
We assume the following two scenarios for the colli- 
sions: (i)In+In — ► quark gluon plasma (QGP) — > 
mixed phase of quarks and hadrons — > hadronic 
phase and (ii) In+In — * hadronic phase and check 
by comparing with data which is the possible sce- 
nario realized in these collisions. In contrast to ear- 
lier studies [1, [t| [HI H3 we focus on the sensitivities 
of the slopes on the equations of state and transition 
temperature, T c . 

The rate of thermal dilepton production per unit 
space-time volume per unit four momentum volume 
is given by[l Hi 
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where a is the electromagnetic coupling, ImlK is 
the imaginary part of the retarded photon self en- 
ergy and /be (E, T) is the thermal phase space factor 
for Bosons, where the energy E should be replaced 
by u^Pn for an expanding system having four ve- 
locity m m at each space-time point of the system, 
p 2 (= Pfj.p 1 ' 1 ) = M 2 is the invariant mass square of 
the lepton pair and 
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arises from the final state leptonic current involving 
Dirac spinors, m in Eq. [2] is the muon mass. 
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The major sources of thermal dileptons are quark- 
anti-quark annihilation in the QGP phase and the 
decays of light vector mesons (p, to and cj>) in the 
hadronic phase. The other non-thermal sources are: 
(i) Drell-Yan process originating from the interac- 
tions of quarks and anti-quarks of the colliding nu- 
clei and the decays of various mesons (e.g. n,T], 
u>, p, T)', 4> etc) after the fireball freeze-out. As 
the non-thermal contributions have been subtracted 
from the data under consideration we concentrate 
only on the thermal emissions. To evaluate the 
thermal dilepton spectra the static rate of emission 
(dR/ d 2 m,TdM 2 dy) has to be convoluted with the 
space-time dynamics governed by the relativistic hy- 
drodynamics as follows: 
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d 4 x is the 4 dimensional volume element, tot is the 
transverse mass and y is the rapidity. The limits 
for integration over invariant mass, M can be fixed 
according to the experimental measurements. For 
the experimental data for M m ; n < M < M max , the 
transverse mass mx is defined as mj- = y/ M 2 + p 2 ^, 
where M = (A/ min + (M max )/2. 

The emission rate from the QGP has been ob- 
tained from the following processes qq — > l + l~, 
qq -> gl+l~ gq{q) -> q{q)l + l~) (see Ref. [H, 
for details). 

The following parametrization [f| [l4[ has been 
used to evaluate the dilepton emission rates from 
light vector mesons (p, lo and <fi): 
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The parameterizations of the vector mesons spectral 
function are consistent with the experimental data 
from e+ e~ — > V(p, uj or 4>) processes [!, H, [3 . fv 
is the coupling between the electromagnetic current 
and vector meson fields, my and Ty are the masses 
and widths of the vector mesons and ujq is the contin- 
uum threshold above which the asymptotic freedom 
is restored. We have taken a s = 0.3, S — 0.2GeV, 
luq = 1.3 GeV for p and to. For <p we have taken 
luq = 1.5 GeV and 8 = 1.5 Gev. The electromag- 
netic current in terms of p, lu and <f> field can be 
expressed as = JP + J^/3 - J<f/3. Therefore, 



the contributions from oj and </> will be down by a 
factor of 9 compared to p. As we have included the 
continuum in the vector mesons spectral functions 
four pion annihilation process [l5[ is not considered 
here to avoid over counting. 

The space time evolution of the system has been 
studied using ideal relativistic hydrodynamics with 
longitudinal boost invariance [161 ] and cylindrical 
symmetry (l7j . The initial temper ature(Ti) and 
thermalization time (rj) are constrained by the fol- 
lowing equation [l8| for an isentropic expansion: 
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where, dN/ dy= hadron multiplicity, Ra is the effec- 
tive radius (evaluated by using the formula Ra ~ 

1/3 

Npart) 01 the system, £(3) is the Riemann zeta 
function and a = ir 2 g/9Q (g = 37 in the QGP), 
A^=number of flavours. The initial radial velocity, 
Vr(i~i, r) and energy density, e(ji, r) profiles are taken 
as: 



v(Ti,r) = 
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where the surface thickness, 6 — 0.5 fm. So far 
there is no consensus on the value of T c , it varies 
from 151 MeV [H to 192 MeV Q. In the present 
work we assume T c — 175 MeV. Results for a higher 
T c will also be discussed. In a quark gluon plasma 
to hadronic transition scenario - we use the bag 
model EOS for the QGP phase and for the hadronic 
phase all the resonances with mass < 2.5 GeV have 
been considered [2l| . The transition region has been 
parametrized as follows 23]: 

s = f(T)s q + (1 - f(T))s h (8) 

where s q (s^) is the entropy density of the quark 
(hadronic) phase at T c and 

/(D = i(l+tanh(^^)) (9) 

the value of the parameter T can be varied to make 
the transition strong first order or continuous. 

The ratios of various hadrons measured experi- 
mentally at different ^/snn indicate that the the sys- 
tem formed in heavy ion collisions chemically decou- 
ple at a temperature (T c h) which is higher than the 
temperature for kinetic freeze-out (Tf) determined 
by the transverse spectra of hadrons [2J]. There- 
fore, the system remains out of chemical equilibrium 
from T c h to Tf. The chemical non-equilibration af- 
fects the dilepton yields at two levels: a) the emis- 
sion rate through the phase space factor and b) the 
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FIG. 4: Same as Fig. H for 0.6 < M (GeV)< 0.9 



space-time evolution of the matter through the equa- 
tion state. The value of the chemical potential and 
the its inclusion in the EOS has been taken in to 
account following [25[ . It is expected that the chem- 
ical potentials do not change much for the inclusion 
of resonances above A. 

To evaluate the tut distribution of dilepton origi- 
nating from the decays of vector mesons after freeze- 
out the Cooper-Frye formula [26| has been used. For 
a special case of unstable vector mesons we need to 
know the thermal phase space factor corresponding 
to an unstable Boson which is given by 



,/unstablc 
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where po = \/p 2 + M 2 , g is the statistical degen- 
eracy and p(M) is the spectral function of the vec- 
tor meson under consideration. For stable particle 
p(M) reduces to a Dirac delta function and conse- 
quently the usual phase space factor for a stable par- 
ticle is recovered. Therefore, the tot distribution of 
dimuons from vector meson decay after the freeze- 
out is given by 
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FIG. 5: Same as Fig. H for 1.0 < M (GeV)< 1.4. 
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where Fy* is the total decay width of the vector 
meson, V. The contributions from Eqs . l3l and fTTI are 
added for the description of the mr spectra of the 
data. 

Before discussing the tut spectra it is important 
to mention here that initially there was no collective 
flow (v(ri,r) = 0), the entire energy of the system 
was thermal. With the progress of time some part 
of the thermal energy get converted to the collec- 
tive (flow) energy for a system undergoing hydro- 
dynamic expansion. The measured mj spectra of 
muon pairs therefore, contains contributions from 
both thermal as well as collective degrees of freedom 
i.e. the inverse slope parameter, T e g can be written 
as T e ff = T + My 2 , as mentioned earlier. Therefore, 
it is important to know the domains of M where 
the thermal contributions from early (quark matter) 
and late (hadronic) phases dominate corresponding 
to small and large radial flow respectively. 

In Fig.l, the invariant mass distribution of dilep- 
tons has been shown for a temperature T = 0.2 GeV. 
It is clear that the quark matter dominates over the 
hadronic matter contributions in the mass windows 
below the p-peak and above the (f> - peak. It is ex- 
pected that the slope parameters extracted from the 
transverse mass distribution of lepton pairs for mass 
region below the p and above the 0-peaks should 
be similar to quark matter phase. Therefore, slopes 



at these M region will correspond to the early time 
when the radial follow is small. On the other hand 
the contributions in the region of p mass is over- 
whelmingly from the hadronic phase and hence the 
slope at this region correspond to the late time con- 
taining larger radial flow. 

First we discuss results for scenario (i) - here we 
have assumed that a thermalized state of quarks and 
gluons is formed after the collisions which reverts 
to hadronic phase through a very weak first order 
phase transition. The value of T in eq. [9] is taken 
as 25 MeV. The values of initial temperature Tj = 
200 MeV, thermalization time t.; = 1 fm and the 
transition temperature T c — 175 MeV have been 
assumed. We take the freeze-out temperature Tf = 
135 MeV which can reproduce the slope of the 4> 
spectra measured by NA60 collaboration for In+In 
collisions 27J. With these inputs transverse mass 
spectra of dileptons have been evaluated for different 
invariant mass windows. 

The resulting mj- — M spectra of /x + p~ are com- 
pared with the experimental data obtained by the 
NA60 collaboration for In-In collision for ^/snn = 
17.3 GeV for different mass window in Figs.[2][5j The 
present calculation agrees well with the data which 
is shown as solid curves for all the mass ranges. 

Finally in Fig[5]the effective temperature obtained 
from the inverse slope of these spectra have been 
plotted and compared with the data. The slopes 
have been estimated from theoretical results (shown 
by solid lines in Figs. [5]© by parameterizing to an 
exponential function within the (mj- — M) range 
0.3 < m T - M(GeV)< 1.0. It is clear from the re- 
sults that the slope at high M region is reproduced 
well if the source at high M region is predominantly 
partonic. A similar non-monotonic behaviour is ob- 
served in the variation of the elliptic flow fa) of 
photons as a of transverse momentum [28l . l29j | . 

Next we discuss results for scenario (ii). In this 
scenario we evolve a hadronic matter of initial tem- 
perature, Tj=200 MeV up to freeze-out temperature 
Tf = 135 MeV with an equation of state which in- 
cludes mesons and baryons up to mass 2.5 GeV. The 
equation of state obtained for such a scenario is soft 
compared to an initial QGP state, the average ve- 
locity of sound (c s ) is c 2 s — 0.16. With this scenario 
we fail to reproduce the data. 

Keeping the values of % = 200 MeV and T f = 
135 MeV fixed as above we increase the value of T c 
to 192 MeV and evaluated the slopes again for all 
the invariant mass windows shown in Figs. [2EJ In 
this case the contributions from the QGP phase is 
considerably reduced so the slopes of the tot spectra 
is pre-dominantly determined by the hadronic phase 
for all the invariant mass windows. The slopes for M 
beyond the <fi peak is not determined by the (early) 
QGP phase but by the (late) hadronic phase. As a 
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FIG. 6: The effective temperature obtained from the 
inverse siope of the invariant mass spectra. The solid 
line indicates results for T* = 200 MeV, T c = 175 MeV 
and Tf = 135 MeV. 



result we observe a monotonically increasing trend 
of T e ff as a function of M. This also indicates that 
the source of lepton pairs at large (M ~ 1.2 GeV) 
is partonic. The sensitivity of T c g to T c for various 
M windows may be used to constrain the transition 
temperature. 

In summary we find that the experimentally mea- 
sured slope of the dilepton spectra by the NA60 
collaboration can be explained theoretically if one 
assumes the formation of a partonic phase initially 
which reverts to hadronic phase through a weak first 
order or continuous transition. A strong first order 
transition does not reproduce the data. A hadronic 
initial state with realistic equation of state fails to 



reproduce the data. It is further observed that the 
spectra is sensitive to the value of the transition 
temperature. A larger value of T c (closer to Ti in 
this case) fails to reproduce the non-monotonic be- 
haviour of the slope, because in this case most of 
the contributions for all the invariant mass windows 
to the lepton spectra come from the late hadronic 
phase resulting in monotonic increase of T e g with 
M. Therefore, the study of T e g may be used to es- 
timate T c . 

Currently the modification of hadronic spectral 
functions in a thermal bath and its effects on 
electromagnetic radiation is a field of great inter- 
est. The invariant mass distribution of lepton pairs 
are sensitive to both the pole shift and broaden- 
ing @, i, H M Hi, H, H S M p. But 

the pt spectra of the EM radiation is insensitive 
to the broadening of the spectral function provided 
the integration over the M is performed over the 
entire region. This is because broadening does not 
change the density of vector mesons significantly (see 
also [33(| ) . However, the number density of vector 
mesons depends on the nature (shape) of the spec- 
tral function within the integration limit. Therefore, 
the pt spectra may change due to broadening when 
the integration over M is done in a limited M do- 
main. By introducing the broadening of p through 
a toy model T p = T° p (l - T 2 /^)- 1 (where T° = 150 
MeV is the width of the p in vacuum) we have ob- 
served that the resulting change in T c q is small. 
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